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We construct simple analytic models of the 5-matrix, accounting for both scattering resonances 
and smooth background contributions for collisions that occur below the s-wave threshold. Such 
models are important for studying confinement-induced resonances such as those occurring in cold 
collisions of ^^^Cs atoms in separated sites of a polarization-gradient optical lattice. Because these 
resonances occur at negative energy with respect to the s-wave threshold, they cannot be studied 
easily using direct numerical solutions of the Schrodinger equation. Using our analytic model, we 
extend previous studies of negative-energy scattering to the multichannel case, accounting for the 
interplay of Feshbach resonances, large background scattering lengths, and inelastic processes. 

PACS numbers: 11.55.Bq,37.10. Jk,34.50.-s,02.30.Mv, 



I. INTRODUCTION 

The ability to control ultra-cold atom-atom interac- 
tions has opened the door to a wide variety of funda- 
mental and applied studies, including the production of 
ultra-cold molecules Juja, [3, S] , simulations of condensed 
matter phenomena [5|, |6| , and quantum information pro- 
cessing '7]. The tools that have been central to this de- 
velopment include designer atom traps, for example, op- 
tical lattices and controllable scattering resonances 
such as a magnetic Feshbach resonance @ . Both of these 
can be used to manipulate the two-body scattering pro- 
cess, thus affecting the strength of the interaction, the 
nature of the resulting two-body states, and more general 
many-body phenomena. Examples include confinement- 
induced resonances p^ . bound-states with repulsive in- 
teractions and Feshbach resonances in band struc- 
tures [111. 

A particular example that we have explored previously 
is a trap-induced resonance (TIR) that occurs as a re- 
sult of interaction between atoms that are confined to 
spatially separated harmonic traps [H, This hap- 

pens as a consequence of a molecular bound state that 
becomes resonant with the vibrational state of the sepa- 
rated atoms due to a quadratic rise in the light shift when 
the two atoms approach one another, as shown schemat- 
ically in Fig. [T] A strong resonance can occur when the 
confinement of the wave packet in the trap is on the or- 
der of (or smaller than) the free-space scattering length. 
This s-wave resonance is analogous to a higher-partial 
wave shape resonance occurring in free space, but here 
the tunneling barrier arises from the trap rather than 
from an angular momentum centrifugal barrier. Because 
of this tunneling, the interaction occurs at "negative en- 
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ergy" values with respect to the free-particle s-wave scat- 
tering threshold. 
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FIG. 1: Schematic of the effective potential between atoms 
trapped in separated wells of an optical lattice, as a function 
of the relative coordinate in the direction of trap separation 
(not to scale). At short range there is a molecular binding 
potential. At long range the relative coordinate is bound by 
the traps. For a given separation and well depth, one of the 
trap vibrational states may become resonant with a bound 
state of the two-body interaction potential, resulting in a trap- 
induced resonance. 

In previous work on the TIR in cesium , whose large 
scattering length gives rise to a strong resonance. Stock 
et al. extracted the scattering length at negative energy 
for the single channel case by an explicit integration of 
the radial Schrodinger equation at negative energies us- 
ing the Numerov method. Such a procedure has limited 
utility; the solutions are unstable since the wave func- 
tion blows up in the tunneling barrier. The situation 
gets very complicated as soon as there is more than one 
channel. A proper numerical technique must ensure that 
open channels are propagated along with the exponen- 
tially decaying closed channels while maintaining accu- 
racy to relevant digits. Most coupled channel codes that 
incorporate such situations (i.e. propagating closed and 
open channels), eventually drop the closed channels be- 



2 



h i j k 1 m n o p 
mp _4 _3 _2 -1 1 2 3 4 



F=4 



F=3 rrr 

mp -3-2-10123 
g fed c b a 

FIG. 2: Energy level diagram of the 6S1/2 hyperfine manifold 
of ^''■^Cs (not to scale). The solid (dashed) lines correspond to 
the levels in the absence(presence) of the an external magnetic 
field. The magnetic quantum numbers and the corresponding 
level-labellings are shown. 



yond certain radius since they are only interested in open 
channels. On the other hand, the problem of TIR consid- 
ered in this paper requires us to integrate to large enough 
radius for determining a good asymptotic log derivative 
(for both open and closed channels). We emphasize that 
we are not allowed to drop any channels, since in the end, 
we are required to extract the asymptotic log derivatives 
for both open as well as closed channels. 

To remedy this, we consider here analytic models of 
the multichannel S'-matrix. Simplified analytic models 
have been employed in previous studies of ultra-cold col- 
lisions and scattering resonances. Julienne and Gao have 
predicted the properties of Feshbach resonances based on 
the analy tic p roperties of the van der Waals long range 
potential [l5|. Marcelis et al. have used analytic mod- 
els to describe the interplay of open and closed channels 
in the context of Feshbach resonances associated with 
a large background scattering length [l^. These ana- 
lytic models, while simple in nature, are able to encap- 
sulate the necessary physics in just a few parameters. 
These parameters can then be incorporated into build- 
ing model many-body Hamiltonians, an example being 
the two-channel model used for describing resonance su- 
perfluidity in a two component Fermi gas [iTj . 

Our goal in this article is to develop an analytic model 
that can be used to study the TIR for Cs atoms trapped 
in an optical lattice. In Sec. II we review the basic 
physics that gives rise to the TIR, its application in Cs, 
and show the limitations of direct numerical solutions, 
even for single channel scattering. Section III contains 
the heart of our new results. We review the basic reso- 
nant scattering phenomena and how they are modeled 
analytically in the S'-matrix. We then apply this to 
determine expressions for the negative-energy scattering 
length in a nontrivial multichannel scattering process, 
relevant to an experimental observation of the TIR. We 
summarize our results in Sec. IV. 



II. SCATTERING RESONANCES IN "^CS 

We consider the scattering of two ^'^^Cs atoms in their 
6S1/2 electronic ground state, trapped in an optical lat- 
tice. The Zeeman hyperfine structure of this manifold 
is shown in the Fig. [2] with magnetic sub-levels labeled 
as a,b,...p for convenience. Henceforth, all two-atom 
scattering channels will be denoted by the relevant pair 
of these sub-levels. To begin with, we consider the scat- 
tering in the \ap) channel. This is motivated by studies 
of controlled collisions via spin-dependent transport in 
polarization gradient lattices [l^. In these spin states, 
two atoms that are separated by A/4 in a lin-perp-lin op- 
tical lattice can be transported into the same well in a 
lin-parallcl-lin optical lattice via a rotation of the laser 
polarization. By angular momentum conservation, be- 
cause these are "stretched states", and ignoring small 
spin-dipolar and second-order spin-orbit coupling p^ . 
the dominant s-wave collision does not couple this chan- 
nel to any other channel. The result is an elastic phase 
shift that can be used to implement an entangling two- 
qubit logic gate. In addition, as the wells approach one 
another, there will be a TIR that can strongly affect the 
two-atom interaction p^ . 

The properties of the TIR follows from a simple model 
of the two-atom system. We express the Hamiltonian for 
the atoms in the \ap) scattering channel in center-of-mass 
coordinates, R, and relative coordinates, r, as 

HcM ^^ + Im.^R^ (1) 

Hroi = ^ + l^ILo^\r-A2\^ + V{r), (2) 

where fi is the reduced mass, Az is the separation of the 
traps, and V{r) is the inter-atomic potential. In prin- 
ciple, the TIR can be seen by diagonalizing the above 
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FIG. 3: Scattering length for the \ap} channel in units of Bohr 
radius a_B as a function of scattering energy E. 
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Schrodinger equation using the precise Cs2 molecular po- 
tential projected on \ap) for V{r). This is a non-trivial 
task, however, since there is a huge separation of length 
scales between the molecular potential and the exter- 
nal trapping potential, and the displacement of the trap 
from the zero of the relative-coordinate makes the system 
anisotropic. Instead, we treat the molecular potential 
through a contact pseudo-potential 20], 

V{r,E) = '^-^a{E)5{T)^. (3) 
/i or 

Here a{E) is the energy-dependent s-wave scattering 
length, determined by direct numerical integration of 
the Schrodinger equation based on the s-wave scattering 
phase shift of the known Cs2 molecular potential in the 
absence of the trap, according to a{E) = — ta,n5o{E)/k. 
The energy is then chosen self-consistently to solve the 
Schrodinger equation, including both the boundary con- 
ditions at short-range due to the atomic interaction and 
at long-range due to the trap [l^l- For sufficient sepa- 
ration between the traps, the lowest energy eigenstates 
drops below the threshold of the molecular potential. 
Thus the "negative energy" scattering states that are in- 
accessible in free space, become opened by the trapping 
potential. 

The scattering length at positive energies for the \ap) 
channel, calculated using a numerical solution to the ra- 
dial Schrodinger equation based on the well-established 
^^•^Cs dimer potential, is shown in Fig. [3l A resonance 
exists Sui E = 4.03 jiK (here and throughout, energy is 
measured in temperature units) due to a bound state 
very close to zero energy. Finding the scattering length 
at negative energies via equivalent numerical integration 
is highly unstable, as the wave function blows up in the 
classically forbidden region. Even if one manages to do it, 
it is necessary that the integration be sufficiently stable 
for large r, well into the asymptotic region of the poten- 
tial, in order to extract a meaningful scattering length 
[l^ . In addition, because we are in the neighborhood 
of a scattering resonance near zero energy, the strong 
variation of the scattering length with energy makes the 
requirement for a robust numerical solution even more 
demanding. To address these problems, we develop ana- 
lytic models that will allow us to calculate the scattering 
matrix below threshold. In doing so, we will extend the 
method to a more complicated process that occurs for 
the multichannel scattering case. 

III. ANALYTIC MODEL 

Near-threshold scattering is dominated by resonant 
phenomena. Such resonances can arise from a variety 
of different physical mechanisms, and thereby affect the 
form of the analytic model. We identify the nature of 
the resonance based on our understanding of the phys- 
ical processes and the location of the poles in the 5*- 
matrix. Away from resonance, the scattering properties 
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FIG. 4: Energy dependent phase shift 5{E) for s-wave scat- 
tering on the \ap) channel. The circles represent numerical 
data from coupled channel calculation, the green curve is the 
analytical fit using Eq. (0 , and the red curve is the analytical 
fit using Eq. p. 

are smooth functions and therefore can be modeled by 
a few free parameters. The total S'-matrix thus factors 
into resonant and background contributions, 

S{k) = 5bg(fc)5rcs(fc). (4) 

For the scattering on the single channel \ap), the res- 
onant behaviour can arise only from a bound state or 
a virtual-bound state near threshold; any other type of 
resonance such as a Feshbach will require the inclusion 
of other channels. A direct numerical integration gives 
a bound state with a binding energy Ei, = 246.2 riK . 
Thus, the resonant part of the single channel S'-matrix 
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FIG. 5: Analytic continuation of the scattering length for 
the \ap) to negative energies. The circles represent data from 
numerical integration of the Schrodinger equation. 
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can be written as [2lj, 



IK^ 



k 



(5) 



where = Eb, representing a pole of the ^-matrix 

in the fc-plane on the positive imaginary axis. There is 
no need to single out the other bound states, as they are 
energetically too far away, and their effect is absorbed in 
the background part. This remaining part can be written 
in the low energy limit fc — )■ as, 



Shgik) — exp(— 2iabgfc), 



(6) 



where Obg is the background scattering length that en- 
capsulates the effect of all other non-resonant processes, 
including other deeply bound states. Now we can write 
the complete S'-matrix element analytically with just one 
free parameter, Obg, the value of which can be determined 
by fitting one positive low energy point to the equation 



(5(fc) = — flbgfc — tan ^{k/Kb). 



(7) 



In Fig. m we plot the scattering phase shift as a function of 
the scattering energy, obtained via the full coupled chan- 
nels calculation and compare it with the one obtained 
from Eq. ([7]). We see good agreement at low energy but 
for energies beyond 10 a slight difference is noticed. 
This is expected since the linear form of the background 
phase shift is only valid at low energies. To remedy this 
problem, we use a higher order expansion for the back- 
ground part given by 



5bg(i?) = 



-l/obg + rok'^/2 + ik 
-l/ahg + rok^/2-ik' 



(8) 



where we have added an additional free parameter, rg, 
the effective range. As before, we determine both the 
parameters by fitting two low energy data points. The 
inset of Fig. [4] shows excellent agreement of this improved 
model with the numerical solution. 

Given the form of the S-matrix, we can predict the 
scattering properties in the \ap) channel at negative en- 
ergy values by performing an analytic continuation of the 
S'-matrix to the imaginary fc-axis. From this one can con- 
sistently define the scattering length at negative energies 
E = -h'^K'^/2n by 



a{iK) 



tan(5(iK)) 



(9) 



In Fig. [Hlwe plot the scattering lengths obtained from the 
analytic procedure discussed above. These agree with the 
direct numerical integration just below threshold [l4|. 

We now turn to a more complex situation: collisions 
between \a) — \F — 3, nip = 3) and |o) = |F = 4, mj? = 
3) or the \ao) channel. This is motivated by the follow- 
ing experimental considerations. Controlled collisions via 
spin-dependent transport in polarization gradient optical 
lattices [l8| is hampered by inhomogeneous broadening 





FIG. 6: Schematic of the potential energy curves for different 
channels. The figure on the left shows potential energy as a 
function of the internuclear separation. The right figure shows 
the scattering threshold for different channels as a function of 
the external magnetic field B. At B = Bq, the bound state 
corresponding to the \bp) channel crosses the \ao) channel 
threshold. The \bp) channel becomes closed for _B > 0. 



arising from unwanted real or fictitious magnetic fields 
(due to elliptically polarized light at the atomic posi- 
tion) [131 . This is a particularly deleterious effect for the 
\ap) states that see a strongly varying difference in their 
optical potentials along the transport. In a lin- angle- lin 
optical lattice at very large detunings, atoms in the \a) 
and |o) states experience almost the same shift due to the 
fictitious magnetic field. Any residual broadening is due 
to the finite detuning effects (giving a differential scalar 
light shift) and the true magnetic field inhomogeneity. 
Controlled collisions in the \ao) channel thus offer the 
advantage of higher degrees of coherence, with potential 
applications to quantum information processing. 

To treat scattering with the incoming \ao) channel, we 
must account for the exchange interaction, which leads to 
spin-changing collisions that preserve the total projection 
of angular momentum along a quantization axis. In this 
case, the s-wave collisions couple the \ao) channel to the 
\bp), |aa), |oo), and \pn) channels. At low energies, small 
compared to the hyperfine splitting, |oo) and \pn) are 
energetically closed. Moreover, in the presence of any 
positive magnetic field B > 0, the channel \bp) shifts to a 
higher energy compared to the \ao) channel, as depicted 
in Fig. [6] For energies that are smaller than this shift, 
\bp) channel is also closed. The movement of this channel 
from below to above the threshold can lead to a Feshbach 
resonance that strongly affects the scattering process, as 
discussed below. 

The scattering phase shift for the \ao) channel is cal- 
culated by a full coupled-channels calculation. In Fig. [7| 
we plot sin^{S{E, B)) as a function of scattering energy 
E and the external magnetic field B in the range of a few 
hundred mC. As is seen in this figure, there is a scatter- 
ing resonance along the dashed line where S{E, B) = tt/2. 
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Also, since the resonance moves monotonically upwards 
in energy as a function of the B field, it is clear that 
this resonance is a Feshbach resonance. We confirm this 
by calculating the bound state energy and find that it 
changes sign at approximately B = 30 mG as shown 
in Fig. [H There also exists a threshold Et{B) corre- 
sponding to the opening of the \bp) channel (shown by 
the white solid line in Fig. [7]). Across this threshold the 
number of open channels changes by one, as reflected in 
the abrupt change in the {E, B) dependence of the scat- 
tering properties. This is indicated in Fig. [5] where we 
plot the scattering phase shift as a function of energy for 
various values of the B field. 

Based on the above understanding of the mechanisms 
that lead to scattering resonances, our goal is to build an 
analytic model that will allow us to predict the scattering 
lengths at negative energies where the TIR is predicted 
to occur. To begin with, we will again assume that the 
diagonal element of the S matrix can be modeled as 

Sao - {ao\S{B, k)\ao) - 5bg(S, fc)5fesh(S, k) (10) 
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FIG. 8: Location of a molecular bound state as a function of 
magnetic field. Data points shown with circles correspond to 
values obtained from direct diagonalization of the two-atom 
Schrodinger equation in a large quantization volume. The 
solid line indicates the position of the bare \hp) channel bound 
state which moves up with B-field. Inset shows the shift of 
the bound state due to dressing by the \ao) channel 



where S'bg is a smooth function describing the back- 
ground contribution and S'fosh is the contribution arising 
from the Feshbach resonance. As before, the resonance 
leads to a pole in the S'-matrix. In fact it can be shown 
that every element of the multichannel S'-matrix has a 
pole corresponding to this resonance [l^. Also, since 
the location of the resonance moves upwards in energy 
almost linearly as a function of the B field, it is fair to 
assume that the bound state solely resides on the \bp) 
channel and is shifted by an amount Afesh due to cou- 
pling to the \ao) channel. Therefore without considering 
the off-diagonal elements of the S'-matrix, the resonant 




FIG. 7: Scattering phase for the \ao) channel as a func- 
tion of energy and magnetic field, plotted as surface plot of 
s\v?{5(E,B)). The solid line corresponds to the boundary in 
the (E,B) plane that separates the region where channel \bp) 
is closed and open. The dashed line corresponds to the points 
where 5{E,B) — n/2. The color scheme used is linear with 
dark red corresponding to the value 1 and dark blue to 0. 



part can be modeled by a pure Breit-Wigner pole, 

c (E ^. E-E,{B)-A,,^^-^T{E)/2 
^^-'^^'''''^-i?-i?.(B)-A,esh + ^^(i?)/2' 

so that that Feshbach contribution to the phase shift is 



ffcsh 



{B, E) = ~ tan" 



T{E)/2 



E - Eb{B) - A 



fcsh 



(12) 



Here Eb{B) is the location of the bare bound state of the 
\bp) channel, r{E) is the width of the Feshbach resonance 
and represents the location of the pole of S{E, B) along 
the imaginary axis in the complex E plane. This model 
is restricted to energies below the \bp) threshold for a 
finite magnetic field, Et{B). The discontinuity in the 
scattering phase shift seen in Fig. [5] shows that a different 
model is required as additional scattering channels are 
opened. We leave that problem for future study. 

To specify our model, we need to find the parameters 
that characterize the Feshbach resonance, F and Afosh- 
This is most easily done by fitting the data to the deriva- 
tive of the phase shift 5{E, B) with respect to B, as de- 
termined by Eq. (fT2|) . 



E',{B)V{E)/2 



{E - EbiB) 



Afe.h)2 



(F(i?)/2)2- 

(13) 

The resonance width and the location of the dressed 
bound state can be obtained by fitting a Lorentzian to 
the numerically calculated values as a function of E. In 
Fig. [To] we plot the values obtained from this procedure. 
The fit to the data agrees well with the expected Wigner 
threshold law, r/2 = Cy/E, where we find C = 2.49 
y/fiK. We neglect the small value of the zero intercept. 
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FIG. 9: Scattering phase shift for the \ao) channel as found 
from a full multichannel calculation, as a function of energy 
for different values of magnetic field B. The crossing of the 
\bp) channel threshold is marked by a finite jump in the value 
of dS/dE. 




FIG. 10: Top figure shows that F varies linearly with ^/E. 
This agrees with the Wigner threshold law. As for the energy 
dependence of Afcsh, it can be assumed to be constant and 
approximated by its value at zero energy. 



typically associated with the presence of the inelastic 
component. 

Finally, from Fig. [TO] we see that for E < 0.2 ^K, the 
level shift acquires a constant value of Afcsh(O) ~ —0.79 
jjK. Thus, we model the scattering phase shift by 



5{E) ~ — abgfc — tan 



(14) 



A more complete model at higher energies requires a de- 
termination of the threshold law for Afosh- Marcelis et al. 
[Fq] have shown for a one channel model or more (if losses 
are neglected) that this law can be obtained from the 
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FIG. 11: Comparison between the analytical model and the 
multichannel data for scattering on the \ao) channel. The 
circles represent values of S obtained from multichannel code 
for B = 0.0287, 0.3, and 0.5 Gauss. The solid lines correspond 
to the analytical model developed in this paper. The fit is 
chosen to be made in the regime where the — bp^ channel is 
closed. 



below-threshold behaviour of the Feshbach bound state 
as a function of the B field. While this is true, a di- 
rect application of their model is not possible here since 
we are dealing with a situation where one of the closed 
channel gets opened. As we are mainly interested in near- 
threshold behaviour in order to extract information be- 
low threshold, we will neglect the functional behaviour 
of Afosh at higher energies. 

As in the case of \ap) channel, the background scat- 
tering length acts as a free parameter that is determined 
by fitting the above model to one low-energy data point. 
In Fig. [Tl]we plot the scattering phase shift S{B,E) for 
three different values of the magnetic field B. The plots 
show a very good agreement of the analytical model to 
the numerical multichannel data. 

Given our model, it is simple to obtain the scattering 
properties such as the s-wave scattering length at nega- 
tive energy (below threshold), by analytically continuing 
the above formula to the positive imaginary axis in the 
complex A:-plane. Figure [T^ shows excellent agreement of 
the scattering lengths obtained using our model with the 
numerical data for positive energies. The predicted val- 
ues at negative energies are also shown which agree well 
with the numerical value obtained just above threshold. 

Critical to the use of the TIR for quantum coherent 
coupling between atoms is a favorable ratio of elastic to 
inelastic scattering processes. Near threshold scattering 
at energies E < Et{B) can couple the \ao) channel to 
one other open channel, \aa). Such transitions will lead 
to loss of the atomic pair, as the hyperfine energy will be 
converted to kinetic energies that are much larger than 
the depth of the trap. By unitarity of the S'-matrix, such 
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inelastic processes imply \Sao\ < 1, which can be modeled 
by an imaginary part of the scattering phase shift. For- 
mal theory of multichannel threshold scattering allows us 
to model the energy dependence of this imaginary phase 
shift [U [U [2^ . Nesbet has shown how to separate the 
threshold behaviour arising from smooth background and 
singularities of poles such as those resulting from virtual 
bound states and Feshbach resonances below threshold 
[26j . For a general multichannel situation consisting of 
M open and n closed channels below and above a certain 
threshold, the off-diagonal elements of the S'-matrix near 
threshold can be written as 

c-. . _ pi(<5bg,r-l-i5bg.r)^. ZKpoio«_ , s 

K -\- ^/^polc 

where 7 is some M x n matrix, i e {1,M}, j € 
and Kpoic is the location of the pole on the imaginary 
k axis. In our case, for B > 0, near the \ao) channel 
threshold, this corresponds to one open channel \aa) and 
one closed channel |ao), the relevant S-matrix element 
has the form 

1^12 1' = 1 - \{ao\S\ao)\-' = Aee-^^^-T^-^, (16) 

+ '*polc 

where is some proportionality constant that can be 
treated as a parameter. The imaginary part of the back- 
ground phase shift is given by the Wigner threshold law 
^hg.i — —ahg,ik where abg.i is the imaginary part of the 
background scattering length. 

The values of the parameters of the model are typi- 
cally obtained by fitting Eq. (|16p to the numerical cou- 
pled channels data once the position of the pole, Kpoic, 
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FIG. 12: Figure shows the analytical continuation of the 
scattering length for \ao} channel to negative energies for 
B — 0.0287 Gauss. There is very good agreement between 
values obtained using the coupled channel code (shown by 
diamonds) and the analytic model (solid line). 
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FIG. 13: Inelastic transitions, 1 — jS'aol'^, plotted as a function 
of energy. A fit of the analytic inelastic model Eq. (|16|) is 
shown by solid line to the coupled channel numerical shown 
by red circles. The dashed line corresponds to the prediction 
for negative energies based on analytic continuation of the 
model. 



is known. Marcelis et al. 16] showed that the exis- 
tence of a virtual bound state leads to a rapid change 
in the shift in the energy of a Feshbach resonance, Afcsh 
in Eq. and derived a simple formula that connects 

the pole location to the shift. From inset of Fig. [H the 
rapid change in Afesh near threshold strongly indicates 
the existence of a such a virtual pole. We cannot, how- 
ever, apply the formula in p^ . since that analysis was 
limited to a two-channel model. Here, we have a more 
complicated interplay of three channels, (|ao), |aa), |6p)), 
since the Feshbach resonance arises from the bound state 
in the \bp) channel. A more detailed analysis is required 
to determine the location of Kpoic from first principles. 
Instead, since sufficient data points are available from 
the coupled channels numerical solution, we treat Kpoie 
as an additional fitting parameter. 

Figure [T51 shows excellent fit of the above model to the 
numerical data. Analytic continuation of this model will 
give a precise determination of the off-diagonal elements 
of the S'-matrix at negative energies. This is shown in 
Fig- [IH for the magnetic field value B — 0.5 G. Since the 
above model for the inelastic component of the S'-matrix 
is known to be valid only in a small energy region close to 
threshold, we have limited ourselves to the energy range 
— O.l/ii^T to Q.lfiK for purpose of illustration. 

As a measure of the ratio of inelastic to elastic pro- 
cesses, we define a general energy-dependent complex 
scattering length, 

Re[a] + iln\[a] = - ta.n[5{E)]/k, (17) 

which can be analytically continued to negative energies. 
In Fig. [131 we plot Im[a]/Re[a]. We see that inelastic 
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FIG. 14: Real and imaginary parts of the scattering length 
for \ao) obtained by analytic continuation of the multichannel 
S-matrix at _B = 0.5 G. Solid lines correspond to numerical 
coupled channels data while dashed line is the analytic con- 
tinuation. 

processes should not dominate, even at negative energies 
sufficiently far from the Feshbach resonance. 

IV. SUMMARY AND CONCLUSION 

In this paper we have discussed the need for an ana- 
lytic model of the ^-matrix in order to extract scattering 
properties at negative energies where numerical methods 
fail. Such negative energy solutions are essential for un- 
derstanding trap-induced resonances that involve atoms 
tunneling into regions of the molecular potential that are 
below the threshold. Beyond improved numerical solu- 
tions, these models give us physical intuition with regards 
to the scattering resonances that are critical to devel- 
oping many-body model Hamiltonians that can help to 
explain ultra-cold atomic phenomena [13, H^j . 

We applied this method to study the collision of two 
"'^^'^Cs atoms in separate harmonic traps, a situation simi- 
lar to that of atoms in separate sites of a polarization gra- 
dient optical lattice. Colliding atoms in spin states that 
are chosen so that they are robust with respect to trap- 
ping inhomogeneities are not necessarily optimal when 
considering the scattering process; they can undergo mul- 
tichannel scattering processes, including inelastic loss. 
We studied a specific example of such scattering - colli- 



sions between |F = 3,M = 3) and |F = 4,M = 3). A 
Feshbach resonance occurs at very small magnetic fields 
(~ 30 mG) due to couphng to a bound state in the 
|F = 3,M = 2)|F = 4,M = 4) channel. To treat this, 
we modeled the S'-matrix element via a smooth back- 
ground component with an imaginary part in its scatter- 
ing length and elastic single-resonance Feshbach model. 
The resulting total scattering length agrees well with di- 
rect numerical multichannel scattering solution at small 
positive energies and extends analytically to negative en- 
ergies well beyond the validity of numerical solutions. 

We have also discussed a model consisting of a few 
parameters in order to describe the off-diagonal element 
of S'-matrix that corresponds to inelastic processes. The 
model agrees extremely well with the coupled channels 
calculation at positive energies above threshold. Ana- 
lytic continuation allows us to calculate the scattering 
length at negative energies below threshold. We note 
that, whereas for free space scattering the reaction rates 
for elastic and inelastic processes just above threshold are 
simply related to the complex scattering length [2^, in 
the negative energy case one must account for the tun- 
neling rate from trap to the molecular potential. We will 
treat this in detail in a future publication. 

Within the framework described in this paper, it is also 
possible to study the threshold behaviour arising due to 
the opening/closing of the \bp) channel. This particu- 
lar threshold is interesting because the channel closing 
can be controlled by use of a very small magnetic field. 
Also since \hp) is the same channel that has the Feshbach 
bound state, opening of this channel leads to the disap- 
pearance of this Feshbach resonance. While Feshbach 
resonances in ultra-cold atomic gases have been thor- 
oughly studied in recent years, our study opens up the 
prospect of studying Feshbach resonances in the vicinity 
of such tunable thresholds and their implications to the 
many-body properties of trapped ultra-cold gases. 
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